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In the experimental context of cold-fermion optical lattices, we discuss the possibilities to approach 
the pseudogap or ordered phases by manipulating the scattering length or the strength of the laser- 
induced lattice potential. Using the Two-Particle Self-Consistent Approach as well as Quantum 
Monte Carlo simulations, we provide isentropic curves for the two- and three-dimensional Hubbard 
models at half-filling. These quantitative results are important for practical attempts to reach the 
ordered antiferromagnetic phase in experiments on optical lattices of two-component fermions. We 
find that adiabatically turning on the interaction in two dimensions to cool the system is not very 
effective. In three dimensions, adiabatic cooling to the antiferromagnetic phase can be achieved in 
such a manner although the cooling efficiency is not as high as initially suggested by Dynamical 
Mean-Field Theory. Adiabatic cooling by turning off the repulsion beginning at strong coupling is 
possible in certain cases. 

PACS numbers: 71.10.Fd, 03.75.Lm, 32.80.Pj, 71.30.+h 
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I. INTRODUCTION 

One of the most exciting possibilities opened by re- 
search on cold atoms in optical traps is to study in a 
controlled manner model Hamiltonians of interest to con- 
densed matter physics. For example, high on the list 
of questions that can in principle be answered by these 
model systems is whether high-temperature supercon- 
ductivity can be explained by the two-dimensional Hub- 
bard model away from half-filling. P, Q As a first step 
towards achieving this goal, the antiferromagnetic phase 
expected at half-filling offers an easier target state that 
occurs at higher temperature in the phase diagram of 
cuprate superconductors. 0, 0] 

In optical lattices, the two spin species occurring in 
the Hubbard model are mimicked by atoms in two differ- 
ent hyperfine states. The cooling of these atomic gases 
necessary to observe ordered states has been discussed 
before, jl, [f| It has been recently pointed out however 
that there is an additional mechanism, [fjj] akin to Pomer- 
anchuck cooling in liquid Helium 3, that is available to 
help in achieving the temperatures where antiferromag- 
netism can be observed. In this mechanism, the temper- 
ature can be lowered by turning on at constant entropy 
what amounts to interactions in the Hubbard model (see 
Ref.(0) for a review). The original calculations for this 
effect were done for the Hubbard model using Dynami- 
cal Mean-Field Theory (DMFT). @ While it is expected 
that this approach will give qualitatively correct results, 
accurate predictions are necessary to achieve the practi- 
cal implementation of this cooling scheme. In the present 
paper, we present such quantitative predictions for the 
isentropic curves of both the two- and three-dimensional 
Hubbard models. We display the results in the usual 



units for the Hubbard model, but also in the conventional 
units used in the context of cold atom physics. 

Solving the problem for both two- and three- 
dimensional lattices fulfills several purposes. First, the 
two-dimensional case is interesting in its own right even 
if long-range order cannot be achieved at finite tem- 
perature in strictly two dimensions (because of the 
Mcrmin-Wagner-Hohenberg theorem). Indeed, high- 
temperature parent antiferromagnetic compounds have 
a strong two-dimensional character. In addition, even 
though long-range order cannot be achieved, there is a 
two-dimensional regime with very strong antiferromag- 
netic fluctuations that is interesting in itself. In this 
regime, a pseudogap appears in the single-particle spec- 
tral weight that is caused at weak coupling by antifer- 
romagnetic fluctuations that have a correlation length 
larger than the single-particle de Broglie wavelength. 

At strong coupling, the pseudogap appears well be- 
fore the long antiferromagnetic correlation lengths occur. 

Also, considering both two- and three dimensions 
sheds additional light on the mechanism for cooling. 

The Hubbard model is defined in second quantization 
by, 



H=- 



(1) 



where c ia (ci a ) are creation and annihilation operators 
for electrons of spin a, n,i a = <J ia Ci a is the density of spin 
a electrons, — t*^ is the hopping amplitude, and U 
is the on-site repulsion obtained from matrix elements 
of the contact interaction between atoms in the basis of 
Wannier states of the optical lattice. [l|, @] We restrict 
ourselves to the case where only nearest-neighbor hop- 
ping t coming from tunneling between potential minima 
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is relevant. In keeping with common practice, t will be 
the energy unit, unless explicitly stated. 

Let us recall the physics of the cooling mechanism pro- 
posed in Ref. [6]. If we denote by / the free-energy 
per lattice site and s the corresponding entropy, then 
s = — (df/dT) u and d = (df/dU) T where d — (n-\ni) 
is the double-occupancy. The chemical potential is kept 
constant in all partial derivatives without further notice. 
In the particle-hole symmetric case at half-filling the den- 
sity is also constant at constant chemical potential. We 
thus have a Maxwell relation 

(£),--(£).■ 

Following Ref. [6(, the shape of the isentropic curves 
s (Ti (U) , U) — est can be deduced by taking a deriva- 
tive of the last equation and using the Maxwell relation 
Eq.© 

where c(Tj) = T {ds/8T) u is the specific heat. If 
{dd/dT) v is negative at small U, then (dTi/dU) s will 
be negative and hence it will be possible to lower the 
temperature at constant entropy by increasing U. Gen- 
erally, double occupancy increases with temperature so 
(dd/8T) u is positive, but it does happen that (dd/8T) u 
is negative, leading to a minimum at some temperature. 
This result may seem counterintuitive. Indeed, at strong 
coupling, namely for interaction strength much larger 
than the bandwidth, such a phenomenon does not occur. 
Double occupancy is already minimum at zero tempera- 
ture. It only increases with increasing temperature. At 
weak coupling however, when the temperature is large 
enough that it allows states to be occupied over a large 
fraction of the whole Brillouin zone, the electrons may 
become more localized than at lower temperature. An al- 
ternate way to understand this minimum is to notice that 
it occurs when the thermal de Broglie wavelength is of the 
order of the lattice spacing: [ll[ At larger temperatures, 
double-occupancy increases because of thermal excita- 
tion while at lower temperature the plane- wave nature of 
the states becomes more apparent and double occupancy 
also increases. The minimum in (dd/dT)u has been ob- 
served in DMFT @,[ll] and also very weakly in Quantum 
Monte Carlo (QMC) simulations of the two-dimensional 
model at U — 4t (see Figure 3 of Ref. [HI]) while in the 
Two-Particle Self- Consistent (TPSC) approach @, ©, [T| 
that we employ along with QMC, very shallow minima 
are observed in three dimensions and are barely observ- 
able in two dimensions depending on the value of U. 
PEE [TtJ As we shall see, in three dimensions a 
minimum in (dd/dT)^ is also predicted by second-order 
perturbation theory. [ll( 

In the next section we discuss the two methods that we 
use, emphasizing the points that are specific to this prob- 
lem. Then we present the results for the constant entropy 



curves in two and three dimensions and conclude. Two 
appendices present calculational details for the entropy 
curves in limiting cases. 

II. METHODOLOGY 

In this section we give methodological details that are 
specific to this work, referring to the literature for more 
detailed explanations of the QMC and TPSC approaches. 

A. Quantum Monte Carlo simulations 

In two dimensions we perform QMC simulations fol- 
lowing the Blankenbecler-Sugar-Scalapino-Hirsch (deter- 
minantal) algorithm. [l8| The standard formula to ob- 
tain the entropy consists in integrating the specific heat. 
However, the evaluation of the latter quantity involves a 
numerical derivative. To avoid differentiating data that 
contains statistical uncertainty, we follow Ref. Q and 
perform an integration by parts to compute the entropy 
from the energy density e 

s(j3,U) = In 4 + /3e (/?,[/) - / e(P',U)d/3' (4) 

J o 

with = l/T in units where the Boltzmann constant 
equals unity. This uses the fact that the entropy at infi- 
nite temperature is known exactly. The integral is calcu- 
lated from the trapezoidal rule on a grid of about twenty 
points spread on a logarithmic scale that extends from 
(3 = to P of order 5 depending on the cases. Each data 
point is obtained by up to 15 x 10 6 measurements for 
the 4x4 lattices and 10 6 measurements for 8x8. By 
comparing with the known result at U = 0, we deduce 
that the error on the integral is of order 2 to 3% at most 
at the lowest temperatures. At large values of U, the 
systematic error due to the discretization of the imagi- 
nary time can be quite large. We checked with U = 14, 
At = 1/10, 1/20 and 1/40 (in units t = 1 which we 
adopt from now on) that At = 1/10 and 1/20 suffice for 
an accurate At — ► extrapolation. 

Size dependence becomes important at low tempera- 
ture. These effects can be estimated from the U = 
case. j!3| The usual formula for the entropy 

a(T,tf = 0) = -^£(/ln/ + (l-/)ln(l-/)) (5) 

k 

with N = L x L (and L even) the number of sites and / 
the Fermi-Dirac distribution, leads to a residual entropy 
at T = given by 

2L — 2 

s(0,0) = — r^ln4 (6) 

which does vanish for L — > oo but which gives important 
contributions for finite L. For example, for L — 4, we 
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FIG. 1: Comparison of TPSC and QMC results for the en- 
tropy as a function of U for different temperatures. Solid 
lines are for TPSC 4x4, triangles for QMC 4x4. The dashed 
line is the TPSC result for the infinite-size lattice limit for 
T = 0.5. Also shown by filled circles for this temperature, are 
the results for QMC 8x8. 



FIG. 2: (Color online) Entropy as a function of 
fiiU /4. From bottom to top, increasing values of U = 
0.5, 1, 2, 3, 4, 6, 8, 10, 12, 13, 14 are displayed. Solid lines are 
for 4x4 lattice and dashed lines for 8x8 lattice. Dotted line 
is the exact atomic limit result for U = 14. 



have s — 0.52, compared with In 4 = 1.386 at T = oo. 
This entropy is easily understood by counting the num- 
ber of ways to populate the states that are right at the 
Fermi surface of the finite lattice in the half-filled Hub- 
bard model. [Ill] At T = 0.3, one can check that the 
relative error between the 4x4 lattice and the infinite 
lattice is about 30% while for the 8x8 lattice it is about 
5%. At T = 0.5, the finite size error for the 8x8 lattice 
is negligible while it is about 5% for the 4x4 lattice. 
Since in this work we concentrate on high-temperature 
results, this will in general not be a problem in QMC. 
The TPSC calculations can be performed in the infinite- 
size limit and for a finite-size lattice. 



B. Two- Particle Self-Consistent Approach 

The TPSC approach has been extensively checked 
against QMC approaches in both two d, and three 
[la ] dimensions. It is accurate from weak to intermedi- 
ate coupling, in other words for U less than about 3/4 of 
the bandwidth, namely U = 6 in d = 2. The double occu- 
pancy is one of the most accurate quantities that can be 
calculated at the first step of the TPSC calculation using 
sum rules. Hence, we can compute the entropy directly 
by integrating the Maxwell relation Eq.© in the ther- 
modynamic limit, or for a finite-size system, using the 
known value of the entropy at U = 0, Eq.lJSJ), to deter- 
mine the integration constant. Earlier results obtained 
with TPSC for the double occupancy may be found for 
example in Refs. [1, [IH, [l6j] . Issues of thermodynamic 
consistency have been discussed in Rcf . 20] . 

Fig. [1] compares the entropy obtained with TPSC 
(solid and dashed lines) and with QMC (symbols) as 



a function of U in d = 2 for different temperatures. 
The QMC calculations are for a 4 x 4 system except for 
T = 0.5 where we also show results for 8x8. The TPSC 
calculations are presented for both 4x4 (solid line) and 
infinite size limit (dashed line for T — 0.5). Down to 
T = 2/3, the results for the 4 x 4 QMC and 4 x 4 TPSC 
agree to better than a few percent for U < 6. One can 
verify from Fig. [T]that at T = 0.5 for U < 6, infinite-size 
limit TPSC and 8x8 QMC results agree remarkably, the 
worse disagreement being less that 10% at U = 6. This 
is expected from the fact that according to the discussion 
of the previous section, finite-size effects are negligible in 
an 8 x 8 lattice in this temperature range. Fig. IIII Bl in 
the following section will compare TPSC estimates of the 
Neel temperature with the latest QMC calculations [2lJ 
in d = 3. There again, U equals 3/4 of the bandwidth 
seems to be the limit of validity. We stress that TPSC is 
in the N = oo universality class |22j so that details may 
differ with the exact result in the critical region. Nev- 
ertheless, it has been checked that even with correlation 
lengths of order 10 or more, the results are still quite 
accurate. 



III. RESULTS 

A. Two dimensions 

The data that is directly extracted from the QMC cal- 
culation is the total energy per site. The entropy ex- 
tracted from this data by numerical integration, Eq.((4]), 
is plotted as a function of (311/4 — U / (4T) on a logarith- 
mic scale for different values of U in Fig. [2] along with 
the exact atomic (single site) limit 



4 




FIG. 3: (Color online) Isentropic curves for d = 2 extracted 
from QMC simulations, including the results of the 8x8 lat- 
tice when they differ from those of the 4x4. Increasing values 
of s are displayed from bottom to top. The first line above 
the horizontal axis is the value of the crossover temperature 
T* determined from TPSC. It stops at strong coupling where 
TPSC ceases to be accurate. 



Satomic (P, U) = In 4 cosh 



pu 



pu 



tanh 



PU 

T 

(7) 

Also plotted in Fig. [2] are the data for an 8 x 8 lattice 
when U = 0.5,1,2, 6. One can check that for p < 1 
(4T > W/2) and U > W, where W = 8 is the band- 
width, the above simple formula describes the data to 
better than 4% accuracy. This is consistent with earlier 
QMC results [13j that found that for U — 10 the specific 
heat above T = 1 is well described by the atomic limit. 
Limiting cases of the above formula are interesting. At 
infinite temperature, /3U = 0, one recovers the expected 
In 4 entropy with the first correction given by 



Satormc (P, U) = 111 4 - {PU) 2 /32 + {{pU)' 



(8) 



In the PU = oo limit, only spin entropy is left so the 
atomic limit result Eq.® reduces to In 2. 

At small values of the entropy, the curves in Fig. [2]at 
small U have a break. This can be understood as a finite- 
size effect given that s ~ 0.5 is the residual entropy for 
a 4 x 4 lattice at U = 0. Lower entropies can be reached 
at larger U without size effects since U lifts the Fermi 
surface degeneracy. [lj| In addition, the results for a 8 x 8 
lattice and small U shown by dashed lines do extend to 
lower values of the entropy. 

Isentropic curves in Fig. [3] are plotted in units of T/t 
and U/t. They are obtained from interpolation of the 
QMC entropy, except for the first line above the horizon- 
tal axis that represents the value of the crossover temper- 
ature T* obtained from TPSC. [23| As discussed above, 
the data in the upper right sector U > 8, T > 1 are quite 
accurately explained by the atomic limit. It should be 
stressed that the slow variation of entropy with T and 



U translates into inaccuracies in the interpolation of the 
isentropic curves that can reach about 10% in this regime. 
When U < 4T and W < 8T, one would expect that the 
high-temperature perturbative result 

s (P, U) = s (P, 0) - i (PU) 2 + O ((/3C/) 4 ) + O (P 3 U 2 W) 

(9) 

derived in the appendices should describe well the QMC 
data. In fact, the term O (P 3 U 2 W) in the range of 
temperatures shown seems to be large enough to essen- 
tially cancel the effect of the leading (PU) 2 term. The 
QMC isentropic curves leave the U — axis with essen- 
tially a zero curvature and are extremely well described 
by the non-interacting result s (P, 0) in Eq.©. More 
specifically, for T > 1, U < 4 the difference between 
QMC and s (/3,0) is less than 3%. At T > 1 again, the 
crossover between the atomic limit value Eq.([7]) and the 
non-interacting value Eq. ([5]) occurs around U = 6 where 
both results differ at T = 1 by about 10% from the QMC 
results. That regime does not lead to an isentropic de- 
crease in temperature concomitant with an increase in U. 
In the non-trivial regime where the entropy may fall with 
increasing U according to DMFT, [|| it is known quite 
accurately that for U = 4, the pseudogap regime where 
antiferromagnetic fluctuations are large begins around 
T = 0.22. Fig. [3] shows that [f| contrary to the three 
dimensional results of the following section, it does not 
appear possible to lower the temperature substantially 
by following an isentropic curve from the U = limit. 
Only a small effect is observed. Even near T = 0.5, 
the isentropic curve deviates only a little bit from the 
non-interacting value but it is quite close to it up to U 
about 6 where a slight downturn in the isentropic curve 
occurs. 24| The flat behavior observed in this regime is 



confirmed by TPSC calculations: the disagreement be- 
tween the two methods down to T = 0.2 is inside the er- 
ror induced by /3-integration of QMC results. The nearly 
horizontal isentropic is not surprising given that the min- 
imum in the temperature dependent double occupancy 
found earlier is shallow in both TPSC [ill, EE El, E| and 
QMC [13, HO, E3 calculations. A very small minimum 
in d (T, U — 4) has been found by extrapolating double- 
occupancy (local moment) to the infinite size limit in the 
QMC calculations of Ref. [l3[ . Note that entry into the 
pseudogap (fluctuating) regime corresponds to a rapid 
fall of d as T decreases. [Ill El El El E3 In fact the 
anticipation of this downfall seems to interfere with the 
formation of the minimum found in higher dimension. In 
the regime where d decreases rapidly as temperature de- 
creases, temperature should increase with U along isen- 
tropic curves, going in the direction opposite to the one 
that would be useful for cooling from the non-interacting 
regime to the fluctuating phase. 

From the large U region, an isentropic decrease in 
U may also lead to a decrease in T, as is obvious al- 
ready from the atomic limit result Eq. ([?])■ In the two- 
dimensional case considered here, the results of the QMC 
calculation in Fig. [3] show that for s < 0.6, adiabatic 
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FIG. 4: (Color online) Isentropic curves for d — 2 (fcs = 1) 
extracted from QMC simulations, expressed in experimental 
units (see text). 



cooling from large U is not possible. All the tempera- 
tures along the s = 0.6 isentropic curves are above the 
fluctuation regime. Hence that regime apparently cannot 
be reached along a single adiabatic curve starting from 
large U and large T. Note however the s = In 2 isentropic 
curve in Fig. [3J It corresponds to the high temperature 
spin entropy in the large U limit (U > T > At 2 /U). [26| 
If one can trap only one of the two atomic species per 
lattice site at random in the large U regime, we are in 
the In 2 entropy case. The lowest temperature that can 
be reached by following this isentropic curve is about a 
factor two above the maximum temperature where the 
strongly fluctuating (pseudogap) regime occurs. 

We now discuss the isentropic curves in terms of ex- 
perimental parameters and units. In optical lattices, 
lasers create a periodic potential defined by a period 
a = A/2, ( A is the laser wavelength) and a depth Vq. 
The energy unit conventionally used in this context is 
the recoil energy: Er = 2ff J| , where m is the mass 
of the fermion. To create a two-dimensional optical lat- 
tice, there is a third standing wave confining the 2D sys- 
tem with a depth that is large enough to prevent out-of- 
plane tunneling. This leads to a 2D on-site energy U re- 
lated to the geometric mean of the confinement strengths 
U/E R = 4V2^(a s /\){V ± /Eg) 1 /' i (Vo/E R ) 1 / 2 where a s is 
s-wave scattering length. [27|, [28| As explained in Refs. 
H, 0, there is a relation to fulfill between a s , Vo and a for 
the one-band Hubbard model to be an accurate descrip- 
tion of cold atoms in optical traps. 

The best way to change only the interaction strength 
for adiabatic cooling is to change the scattering length, 
as can be done by tuning through a Fcshbach reso- 
nance. If only the scattering length is changed, the 
shape of the adiabatic curves will be as in Fig. [3] 
Only the scales need to be changed. All energies in 
that plot are in units (fcs = 1) of hopping t which is 
related to recoil energy and potential strength through 



t = E R {4/^)(Vo/E R ) 3 ^exp [-2^/Vo/E^M 

We can also change U by changing the potential 
strength Vq but clearly this changes also the hopping 
t. Thus for quantitative purposes, we also display the 
preceding isentropic curves in the (Vq/Er,T/Er) plane, 
rather than in the (U/t,T/t) plane. This change in co- 
ordinates is discussed in Appendix [AJ It has a strong 
influence on the shape of the isentropic curves as can be 
seen in Fig. [4] obtained for the value a s /a = 2 10 -3 
and U = E R 4V2n(a s / \)(V ± / Er) 1 / 4 ^/ Er) 1 / 2 with 
V_l/E r = 30. [1| We also display in this figure the 
pseudogap temperature T* determined in the TPSC ap- 
proach. [23| As Vq/Er increases, the system cools down 
along isentropic curves, at least for moderate Vq/Er val- 
ues. For higher values, isentropic curves corresponding 
to s > In 2 eventually bend upwards, while those corre- 
sponding to s < In 2 bend downwards, in such a way that 
there is a large domain of temperatures in the vicinity 
of the large-repulsion spin-entropy value s = In 2. This 
general behavior was also present in Fig. [31 and will be 
seen in 3D too. We checked than the general appear- 
ance does not change for other reasonable values of a s /a 
compatible with the Hubbard model. Given that the 
temperature axis is displayed on a logarithmic scale, it 
thus appears that tuning the potential Vq can be very ef- 
fective in reducing the temperature of the fermions. The 
general cooling trend is due to the decrease in hopping t 
associated to an increase in Vo- Indeed, the ratio of tem- 
perature to bandwidth is constant for isentropic curves 
of non-interacting electrons, hence T decreases monoton- 
ically with decreasing t in this case. This is the mech- 
anism discussed in Ref. @. At Vq < 2.3E R , (not on the 
figure) heating can also occur. @ The presence of inter- 
actions can enhance the cooling compared with the non- 
interacting case. However, cooling down the system 
along an isentropic by increasing Vo does not necessarily 
mean an effective approach of the strongly fluctuating 
regime of the system: indeed, as can be seen from the 
figure, if we increase Vq/Er further than about 12 the 
pseudogap region in experimental units moves away to- 
wards lower temperatures. Note that T* determined by 
TPSC is not reliable at large values of Vq/Er: the limit of 
validity ^ ~ |, corresponds to ~ 12, for our choice 
of a s /a and V±/E R . 



B. Three dimensions 

To discuss the isentropic curves in 3D, let us go back 
for a while to usual units and parameters of the Hubbard 
model. In three dimensions, adiabatic cooling towards 
the antiferromagnetic phase, starting from small U, is 
possible and quite clearly so. In fact, it occurs at high 
enough temperature that perturbation theory (Appendix 
[C| suffices to show the effect. This is made clear by Fig. 
[5] where ds/dU changes sign from negative to positive on 
isothermal curves, as T decreases. By Maxwell's relation 
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FIG. 5: (Color online) Comparisons of full TPSC calculation 
(solid lines) with second-order perturbation theory (dashed 
line) for the entropy as a function of U for different temper- 
atures. 
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FIG. 6: (Color online) Isentropic curves for d = 3 extracted 
from TPSC. Increasing values of s are displayed from bottom 
to top. The lowest solid line is the Neel temperature. Symbols 
are the results of QMC calculations taken from Ref . [211] . 



Eq.©, this reflects the change in sign of dd/dT. When 
the temperature is large enough, the dashed lines from 
second order perturbation theory agree very well, up to 
quite large interaction strength, with the solid lines from 
the full TPSC calculation 

The TPSC results for the isentropic curves in three 
dimensions are exhibited in Fig. IIII Bl In the low- 
temperature regime 2irT/W < 1, TPSC is strictly valid 
only in the U <C W (W = 12) limit. This can be checked 
by comparing the TPSC Neel temperature with that of 
the latest QMC calculations, shown by symbols. [2l[ 
Clearly, the agreement is satisfactory up to U ~ 8 (or 
U ~ 3W/4 as mentioned before) where the Neel temper- 
ature as a function of U saturates according to TPSC 
but begins to decrease according to QMC. 



Despite the fact that TPSC is not valid in the atomic 
limit, it seems to recover the correct result at high- 
temperature even if U > W. Consider for example T = 
1.5. That temperature is reached along the s — 0.8 isen- 
tropic around U = 14 in TPSC and around U = 13 
in DMFT. [f| Similarly, T = 1.5 is reached along the 
s = 0.75 isentropic for U ~ 17 for both TPSC and 
DMFT. The corresponding atomic limit results, that are 
dimension independent, are that s (T — 1.5, U — 11) ~ 
0.80 and s (T = 1.5, U = 14) ~ 0.75. TPSC is closer to 
DMFT than to the high-temperature atomic limit, sug- 
gesting that both approaches take into account the same 
physics at U large in the high-temperature limit. Note 
however that the DMFT and TPSC results are different 
at U = because a model density of states is used in 
DMFT instead of the one following from the exact dis- 
persion relation used in TPSC. 

As in the two-dimensional case, the value of the en- 
tropy at T > 1 and U < W is almost independent of U. 
Contrary to the two dimensional case however, there is 
a region, namely for s < 0.65, where cooling along isen- 
tropic curves down to the interesting regime is possible. 
Cooling however is from T ~ 0.75 to the maximum Neel 
temperature, T ~ 0.4. DMFT predicted cooling to the 
Neel temperature beginning around T ~ 1.1. 

The possibility of adiabatically cooling all the way to 
the Neel temperature starting from large U is also dis- 
cussed in Ref. t 6!]. For this one needs two conditions. 
First, the entropy at the maximum Neel temperature has 
to be larger than the entropy of the Heisenberg antifer- 
romagnet. This condition is satisfied according to Fig. 
IIIIBI since the entropy at the maximum Neel temper- 
ature is around s max = 0.65 while the entropy of the 
Heisenberg antiferromagnet is a constant sh estimated 
in Ref. to be about 50% smaller than In 2. This would 
mean that there is indeed a maximum in the value of the 
entropy at the Neel temperature plotted as a function of 
U. That maximum would be even more pronounced than 
that sketched in Fig. 3 of Ref. [6j. The second con- 
dition to be satisfied is that the temperature should de- 
crease as U decreases along isentropic curves in the range 
sh < s < s max . TPSC cannot tell whether this condition 
is satisfied since for temperatures less than roughly unity 
at strong coupling U > 8 the TPSC results cannot be 
fully trusted. 

Incidentally, if one takes the TPSC results seriously 
up to T ~ 0.7, then it is not possible to cool all the 
way to Neel temperature along the s = In 2 = 0.69 
"infinite" -temperature {U/T <C 1) isentropic curve, con- 
trary to what DMFT suggests, since the minimum in the 
s = 0.7 isentropic curve in Fig. IIII Bl is roughly a factor of 
two above the Neel temperature. As discussed in Ref. [f| , 
DMFT does not give an accurate estimate of the entropy 
at the Neel temperature since the latter is obtained in 
mean-field. (30j This is why the s = 0.7 isentropic curve 
ends at the Neel temperature in that approximation. 

Finally, we discuss experimental units. Once again, 
increasing the scattering length would be the sim- 
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FIG. 7: (Color online) Isentropic curves for d — 3 from TPSC 
calculations, expressed in experimental units (see text). 

plest way to implement directly the Pomeranchuck 
adiabatic cooling discussed in this paper. Indeed 
varying the scattering length value will span the ab- 
scissa axis in Fig lHIBl All energies in that plot 
are in units (k R = 1) of hopping t which is re- 
lated to recoil energy and potential strength through 

t = ^(4/V^)(V£fl) 3/4 exp(-2/ty^).[I3 The 

interaction strength on the other hand is U = 
E R 4V2^(a s /X)(V /E R ) 3 / 4 .^ 

As in the 2D case, we choose to change the potential 
strength Vq, that modifies both hopping t and on-site in- 
teraction U. This is shown in Fig[7]where we display the 
isentropic curves and the TPSC-determined Neel tem- 
perature in the experimental units [311 ]. For the figure, 
we fix a s /a = 2 10~ 3 . Increasing the potential Vq is quite 
efficient to cool down the system in absolute units, but 
the Neel temperature also recedes. Two trends remain: 
for s > In 2 the cooling is not sufficient to reach the an- 
tiferromagnetic region, but a smaller entropy value may 
do. In these units, the isentropic s = In 2 seems "noisy" 
for large Vq/Er values because, in this region, the en- 
tropy surface is flat: this makes the precise location of 
the isentropic curve difficult to determine. 

IV. CONCLUSION 



nism is not very effective, making it impractical to reach 
the low temperature fluctuating regime by this approach. 

In both d — 2 and d = 3 it is possible to adiabatically 
cool starting from the large U regime as suggested in Ref. 
@. However, in d — 2, the fluctuating regime cannot be 
reached along a single adiabatic using this approach. In 
d = 3, although there are encouraging trends, we cannot 
tell unambiguously with TPSC whether the Neel temper- 
ature of the Heisenberg antiferromagnet can be reached 
by decreasing U along a single adiabatic that starts at 
high temperature. 

In the context of cold fermions, changing the strength 
of the laser-induced lattice potential changes both hop- 
ping t and interaction strength U . In units of absolute 
temperature and potential strength then, the shape of 
the adiabatic curves are quite different from those in the 
T/t and U/t units appropriate for the Hubbard model. 
While the change in t can, for some cases, produce dras- 
tic cooling in absolute units, it also changes the shape 
of the lines for the Neel temperature (in 3D) and the 
pseudogap temperature (in 2D). We have plotted the re- 
sults for a given scattering length as examples. Whether 
a given isentropic curve crosses the Neel or the pseudo- 
gap temperature is clearly independent of coordinates. 
To implement directly the type of interaction driven adi- 
abatic cooling discussed in the present paper, one could 
change only the interaction strength by manipulating the 
scattering length with a Feshbach resonance. 

Adiabatic cooling away from half-filling and for other 
cases can be studied with the methods of this paper. 
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TPSC calculations confirm that the physics of adia- 
batic cooling by increase of U (as found in Ref. [6() is 
correct for three dimensions. Our quantitative estimates 
show a smaller but still appreciable effect. Since it oc- 
curs at relatively high temperature, that effect is qual- 
itatively captured already by second order perturbation 
theory. Adiabatic cooling should help in reaching the 
Neel antiferromagnetic transition temperature. Reach- 
ing the antiferromagnetic phase should be a first step in 
the study of d-wave superconductivity in optical traps. 
In two dimensions however, QMC shows that this mecha- 



APPENDIX A: EXPRIMENTAL COORDINATES 
FOR THE TWO-DIMENSIONAL CASE 

The problem of converting from the theoretical units 
(U/t, T/t) to the experimental units (V /E R ,T/E R ) is 
straightforward in three dimensions but it requires some 
discussion in two dimensions. Since [27], [28| 
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and 

Er y/ir \E R/ 
it follows immediately that 




(A2) 




For dcfinitcncss wc choose in this paper as/A = 1CP 3 and 
V±/Er = 30. Clearly, U/t is not a monotonic function of 
Vq- For a given U ft we have two or zero real values of Vq 
as we now discuss. 

For the sake of simplicity, let us define the reduced 
units 



U 

u = — 
t 



1/4 



Then, the relation between u and vq becomes 



-1/2 



exp (2v ) 



(A5) 



(A6) 



(A7) 



This function has a minimum at vq — 1/4. The only 
possible values of u are thus u > 2-y/e and for each such 
values of u, there are two values of Do , one less than 1 /4 
and the other one larger than 1/4. It is the latter that 
we consider as the physical value. Indeed, it corresponds 
to Vq/Er > 1/16 and we know that the Hubbard model 
is valid only for sufficiently large values of Vq/Er. The 
minimum value of U/t for V±/Er = 30 is about 34a s /A 
which is quite a small value. To solve for vo (u) , it suffices 
to rewrite Eq. (|A7[) as 



4 

_ = -4v cxp (-4u ) • 



(A8) 



The solution to y — x exp (x) is the Lambert function Wk 
(also known as "product log" function). Since y (x) is 
non-monotonic, there exists a family Wk of inverse func- 
tions. If U is larger or equal to the bound discussed 
above, then > y > — 1/e and Wk=-i,o can take real 
values. The branch x = W—i (y) has x < — 1, corre- 
sponding to Vq/Er > 1/16 that we want to retain. The 
other branch, x = Wo (y) leads to x > —1 which we do 
not consider here. Hence the solution is 



Er 



-W-i 



u 2 



1 2 



Sir 2 



Er 



1/2N 



(A9) 



APPENDIX B: ENTROPY IN THE LARGE 
TEMPERATURE LIMIT FROM THE 
SELF-ENERGY 

Consider the large Matsubara frequency (equivalently 
high-temperature) limit of the self-energy, [8|] 



S (k,ifc n ) = Uri- 



U 2 n a (1 - n- a ) 



(Bl) 



This formula is valid when ttT is larger than the fre- 
quency range over which T, R " is non-zero. In practice, 
since in all the diagrams that enter the calculation of the 
self-energy, k n = (2n + 1)tvT is compared with £k and 
there is particle-hole symmetry at half-filling, we may ex- 
pect that as soon as the Matsubara frequencies are larger 
than band energies of order ±W/2, namely irT > W/2, 
then the expansion may apply. This is confirmed by the 
numerical results in this paper. The expansion should 
thus be valid for ttT > 4 in d = 2, and ttT > 6 in 
d = 3. If, inspired by the exact atomic result Eq.([7]) we 
replace tt by 4 we recover the limits of validity mentioned 
in the text. In addition to this restriction on temperature 
compared with bandwidth, we note that this asymptotic 
expansion for the self-energy Eq. (|B1[) is clearly a power 
series in U/wT. At half- filling, using the usual canonical 
transformation to the attractive Hubbard model, we can 
see that if we absorb the Hartree-Fock term in the defi- 
nition of the chemical potential then only even powers of 
U enter the expansion, which makes it convergent even 
faster. In addition, it turns out that stopping the ex- 
pansion of £ (k,ik n ) at U 2 /ik n at half-filling reproduces 
the exact result in the atomic limit. Hence, in the spe- 
cial case we are interested in, we expect that this high- 
temperature expansion W/2ttT 1 is excellent for arbi- 
trary values of U, even if strictly speaking it should be 
valid only if U / irT < 1 as well. 

From the above expression for the self-energy and the 
sum-rule, [1] 

§ E E s ( k ' ifc «) G ( k ' ifc «) e " fc "°~ = u ( n T n i> ( B2 ) 

n k 

we can extract the double occupancy d that we need to 
compute the entropy in the high-temperature limit. For 
the Green function, we again assume that W/2nT <C 1. 
This means that we can insert in the previous equation 



G (k,ik n ) = - 



1 



kri 



u 2 



(B3) 



This clearly neglects band effects that would contribute 
to order UW/ (ttT) to double occupancy. The normal- 
ized sum over wave- vectors in the sum- rule Eq. (|B2|) con- 
tributes a factor unity while the discrete Matsubara sum 
can be performed exactly. One finds, 

- ~ tanh (JL\ = U «n T n;) - <n T ) (n L )) . (B4) 
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To extract the entropy, it suffices to use Maxwell's rela- 
tion Eq.@ so that 



s(T,U) = s(T,0) 



g(_i tanh (^)) 



dT 



dU 



s (T, 0) + In 



cosh ( — — 
4T 



U 

—— tanh 
4T 



(B5) 

U 
4T 



The above expression Eq.jB5} with the exact value for 
s(T, 0) neglects terms of order U 2 W/ (irT) 3 . In prac- 
tice, we found that keeping the non-interacting value 
of the entropy s (T, 0) in the above formula does not 



J 



T 



Jr E E S ( k >^«) G ( k < 



-ife„0~ 



n k 



T 



improve the comparison with QMC data in the region 
where W/2irT <C 1 is satisfied, whether U is small or 
large. When we neglect all band effects compared with 
temperature, then s(T, 0) can be replaced by In 4 and 
we recover the atomic limit result Eq.([7|) that can also 
be found from elementary statistical mechanics. It is the 
latter result that is useful to understand the data at large 
values of U. 

Expansion of s (T, U) above in powers of U/4T leads 
to the perturbative result Eq.©. One can also arrive at 
this result by directly neglecting higher powers of U /irT 
and UW/ (irT) 2 in the self-energy and Green function, 



n k 



V 2 n a (1 - 1 



= U{ n] )( ni )-Tj2 



u 2 



-ik n Q~ 



4 (2n + l) 2 (ttT) 2 



U 2 



U (n T ) (nj.) - — = U (n T nj.) 



(B6) 
(B7) 
(B8) 



so that 



s(T,U) = s(T,0) 
- s(T,0) 



u dJrM 

o dT 

U 2 



32T 2 ' 



dU (B9) 
(BIO) 



This result, appearing in Eq.Q, keeps all powers 
in W/2irT, the leading term in U/irT and neglects 
U 2 W/ (7rT) 3 and higher orders (the entropy is an even 
function of U at half-filling) . It does not however assume 
that U/W < 1. It is the large temperature here that con- 
trols the expansion. In practice we found that the above 
formula does not lead to a good description of the QMC 
data in any regime, even small U and large T, unless 
s (T, 0) — > In 4 in the large U regime. This suggests that 

the corrections O are important 

and in fact cancel the leading one. 

Note that in all the results of this section, the dimen- 
sion occurs only in the value of W and in the value of 
s (T, 0) . The atomic limit is independent of dimension. 



APPENDIX C: SECOND ORDER 
PERTURBATION THEORY AND TPSC FOR 
THE ENTROPY 

In the limit U <C W, TPSC reproduces the standard 
perturbative expression for double occupancy. This can 
be demonstrated as follows. In TPSC, double occupancy 



is obtained from the following sum rule and ansatz 0, [l4| 



n - 2(n T n x ) 



U„ 



T 



Xo(q) 



N „ l - hu sp xo(q) 



u 



!«T><«4> 



(CI) 
(C2) 



We used short-hand notation for wave vector and Mat- 
subara frequency q = (q,ig„). Since the self-energy is 
constant in the first step of TPSC, the irreducible sus- 
ceptibility takes its non- interacting Lindhard value Xo(o)- 
In a perturbation theory in U, we can expand the right- 
hand side of the sum rule Eq. (|Cl[) and take U sp = U 
which leads to 

n-2(n T n;) = ^ ^ ( X o(q) + \u^q) 



(n T n;) - (n T )<n;) = ~^U— 22xI(q) 



(C4) 



which shows, as expected, that double-occupancy is 
decreased by repulsive interactions compared with its 
Hartree-Fock value. The above corresponds to the ex- 
pression obtained from direct perturbation theory for 
(nfnj) — (n.f)(n|). The entropy can be obtained, as usual, 
from integration of the Maxwell relation Eq. @ using the 
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known s (T, [7 = 0). This is how the perturbative result 
in Fig. [5] was obtained. 

In the high-temperature limit, we recover results of 
the previous section, as we now proceed to show in two 
different ways. TPSC at the first level of approxima- 



tion obeys the sum rule Eq. (|B2[) that expresses a consis- 
tency between single-particle and two-particle quantities. 
The self-energy in the U/W < 1 and large Matsubara- 
frequency limit has been found in Ref. Q, Eq.(E.lO) 



£ (k,ik n ) = Un_ a + — 



u fu sp + u, 



'ch 



u, 



chn_ a 



(U sp — U c h) 



n1 n i) + 



(C5) 



r 



In the high-temperature limit ttT 3> W/2 that we are 
interested in, the classical (zero-frequency) contribution 
dominates the sum rules used to find U sp and U c h so that 
U sp = U c h = U and one recovers that the l/ik n term has 
the exact form U 2 / (iik n ) used in the previous appendix. 
We thus recover the high-temperature perturbative result 
for the entropy derived there and appearing in Eq. (J9J) . 

Another way to arrive at the same result in a more 
transparent way that uses only the first step of the TPSC 
approach (U/W < 1) is to work directly with the previ- 
ous perturbative result Eq. (|C4[) and evaluate it in the 
high temperature limit. But we first rederive the per- 
turbative result Eq.JUU from Eq.(43) of Ref.0 that is 
valid when the correction of double occupancy from its 
Hartree-Fock value is small, 



(n T n ; ) = (n T ) (n;) 1 + M/ - 



(C6) 



Correcting a factor of 2 misprint in Ref. [8| , the quantity 
A is given by 



A = ilyEE^ (q.<?n) 



(C7) 



iq n q 



with q n a bosonic Matsubara frequency and xo the Lind- 
hard function. Expanding the denominator in Eq. (|C6D 
and substituting n = 1 and U (n-f) (nj.) ~ U/4, we do 
recover the perturbative result Eq. (|C4[) . 

In the limit W/2-kT <C 1, the susceptibility xo scales 
as 1 /q^ which yields terms that are smaller in powers of 
W I (2ttT) than the zero-Matsubara frequency contribu- 
tion. Neglecting these finite Matsubara frequency terms, 



and taking the large W/ (2ttT) limit where / (ek) 
0.5 (1 — 0.5/3ek) , we are left with 



Xo (q,o) 



-E 

q 

I 
2' 



/ (gk) " / (£k+q) 
£k — £k+q 



(C8) 
(C9) 



("T n l> - ("T> ( n l) 1 



^From this at n — 1 we can evaluate that A ~ T j (2T) so 
that, to leading order, the approximate double occupancy 
found from Eq. (|C6|) is 



-)= - - — (cio) 

AT J 4 16T v ' 
which leads again to the high-temperature perturba- 
tive result found at the end of the previous appendix 
Eq. (|B9[) and hence to Eq.@. Even if this time we as- 
sumed U /W < 1 in the last derivation (there is no 
Mott gap in the one-body Green functions), the fact that 
the asymptotic TPSC self-energy Eq. (|C5|) in the high- 
temperature limit reduces to U 2 /ik n plus corrections that 
involve two more powers of U j ttT suggests (but does not 
prove) that the atomic limit is also satisfied by TPSC at 
high temperature. In the high-temperature limit where 
s (T, 0) — * ln4, the result that we just found, Eq.©, does 
reduce to the first two terms of the high-temperature se- 
ries of the atomic limit. A coincidence between atomic 
limit and TPSC was also noted for the attractive Hub- 
bard model in Ref. 1321. 
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